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We study the possibility to create many-particle Schrodinger cat-like states by using a Feshbach 
resonance to reverse the sign of the scattering length of a Bose-Einstein condensate trapped in a 
double-well potential. To address the issue of experimental verification of coherence in the cat-like 
state, we study the revival of the initial condensate state in the presence of environmentally-induced 
decoherence. As a source of decoherence, we consider the interaction between the atoms and the 
electromagnetic vacuum, due to the polarization induced by an incident laser field. We find that 
the resulting decoherence is directly related to the rate at which spontaneously scattered photons 
carry away sufficient information to distinguish between the two atom-distributions which make-up 
the cat state. We show that for a 'perfect' cat-state, a single scattered photon will bring about a 
collapse of the superposition, while a less-than-perfect cat-like state can survive multiple scatterings 
before collapse occurs. In addition, we study the dephasing effect of atom-atom collisions on the 
cat-like states. 

PACS numbers: 03.75.Gg, 32.80.Lg,03.75.Lm 
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I. INTRODUCTION 

Quantum degenerate gases of bosonic and/or fermionic 
atoms have proven very useful tools in the exploration of 
low-temperature many-body quantum physics. In this 
area of research, the initial Bose-Einstein condensate 
(BEC) or degenerate Fermi gas (DFG) serves primar- 
ily as a well-defined starting point for the generation of 
more exotic highly-correlated and/or entangled many- 
body states. The most profound recent experimental 
demonstrations of this approach include the observations 
of the superfluid to Mott-insulator transition in an op- 
tical lattice and the cross-over to a 'fermionized' Tonks- 
Girardeau under quasi-one-dimensional confinement. In 
addition, there has been significant progress towards ob- 
serving the BEC to BCS cross-over in a two-component 
gas of Fermionic atoms P, , and efforts are underway 
to observe the cross-over from a vortex-lattice to a non- 
trivial entangled many-body state in the ground-state of 
a rapidly rotating BEC Hg. 

In the first of the aforementioned experiments, a re- 
versible change from a superfluid to a Mott-insulator 
state was observed when crossing from tunncling- 
dominated to the collision-dominated regime for atoms 
trapped in a lattice potential with repulsive atom-atom 
interactions. The analogous transition for the case of 
attractive interactions involves the crossover from a su- 
perfluid into a state often described as a 'Schrodinger 
cat state' where the atomic population collapses into a 
single lattice site, with the true ground state being a sym- 
metric superposition over all possible lattice sites as the 
final occupied site. In addition to providing fundamental 
insights into the nature of the transition from the quan- 
tum to the classical descriptions of reality, such states 
may have important applications in precision measure- 
ment [B| and quantum information processing. 

In this paper we investigate theoretically the possi- 



bility to use a BEC trapped in a double-well potiential 
to generate highly- entangled Schrodinger cat-type state, 
which can be viewed as a realization of the two-site Bose- 
Hubbard model with attractive interactions. We ad- 
dress important issues such as detection of cat-like states 
and the effects of photon-scattering induced decoherence 
which occurs when far-off-resonant laser-light is used for 
trapping and/or probing the cat states. 

Interest in Schrodinger cat-like states, loosely defined 
as quantum superpositions of macroscopically distin- 
guishable many-body states, goes back to the earliest 
days of quantum mechanics. In 1935, in oder to demon- 
strate the limitations inherent in using quantum mechan- 
ics to describe everyday phenomena, Erwin Schrodinger 
proposed a experiment in the macroscopic state of a 
cat is entangled with that of an unstable nuclei, allow- 
ing the cat to enter a coherent superposition of being 
dead and alive 0. Such states were thought at the 
time to be logically untenable, although today we might 
consider only that they would not be observed due to 
rapid environmentally-induced decoherence @. To this 
day, this famous gedanken experiment illustrates the fun- 
damental difficulties inherent postulating a well-defined 
boundary between a 'classical' level of reality and an un- 
derlying level governed by quantum mechanics, one of 
the fundamental open questions in quantum theory. 

Many of the early gedanken experiments, whose con- 
sideration led to the earliest understanding of quantum 
mechanics, are now being transformed into laboratory re- 
alities. One such effort has been the ongoing quest for 
Schrodinger-cat-like states, with the most striking experi- 
mental results to date being the observation of a superpo- 
sition between clock- wise and counter-clock- wise currents 
in a superconducting quantum interference device and 
the observation of double-slit diffraction for massive Cgo 
molecules . Additional efforts include work in trapped 
ions at low temperature 0,^3, nanoscale magnets 
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and superconductors 0, ITil ITEf . This progress is gener- 
ating deeper insight into the meaning of quantum theory, 
as well as potential applications in the fields of precision 
measurement and quantum information processing. Even 
with the current rapid rate of progress, the quest to ob- 
tain and utilize large cat-like states remains a serious ex- 
perimental challenge, due primarily to the presumed scal- 
ing of fragility (with respect to environmentally-induced 
decoherence) with number of particles. For most applic- 
ations, the utility of the cat-state increases with particle 
number 0, , so there is a strong need to find systems 
in which such maximally-entangled states of large num- 
bers of particles can be generated without rapid collapse 
due to decoherence. 

An atomic BEC [Til IT^ . l2*cj may be such a system, 
as its large degree of spatial and temporal coherence is 
well established by experiment The bimodal BEC, 
which is constructed by isolating a pair of weakly-coupled 
atomic-field modes, is one of the simplest quantum sys- 
tems used in the study of entanglement physics and has 
therefore been the subject of extensive theoretical study 
[22II23I ]. In this system, N bosons are restricted to occupy 
one of two modes, and entanglement is established via 2- 
body interactions between the bosons. For the case of an 
atomic BEC, the two modes can either be two potential- 
wells spatially separated by a potential barrier (so-called 
Double- well BEC's) 0iE3 or two spatially-overlapping 
modes in different hyperfine states |2a , |27| . In addition, 
a recent paper proposes an analogous system in which a 
double- well 'potential' in momentum-space, formed from 
the band-structure in a lattice potential, can be used to 
create cat-states involving two momentum-states [28| . 

The quantum ground-state of the bosonic double-well 
system has been studied for the case of repulsive inter- 
actions in the context of the superfluid to Mott-insulator 
transition (2^, as well as for the case of attractive in- 
teractions. For the latter, it has been shown under cer- 
tain conditions a Schodinger cat-like macroscopic super- 
position state, can be generated |3(J, |3l|. In this cat-like 
many-body state, a large number of atoms are collec- 
tively either in one mode or in another. The dynamics 
of the bimodal system have also been investigated, both 
via a meanfield approach and a full many-body treatment 
within the two-mode approximation. In the mean-field 
treatment, the dynamic evolution is treated as a Bose 
Josephson Junction, which is complementary to that of 
a superconductor tunnel junction [32l l33l| . This treat- 
ment predicts a degenerate ground-state, corresponding 
to all amplitude in one well or the other, but is not ca- 
pable of predicting the dynamical evolution of cat-like 
states, as these states lie outside of the domain of mean- 
field theory. In the two-mode many-body treatment, the 
many-particle cat-like state is predicted to arise during 
the dynamical evolution of an initial minimally-entangled 
condensate state in the presence of repulsive interactions 
[25L l26l | , but with the sign of the tunneling coefficient re- 
versed by imprinting a 7r-phase shift onto one well. In this 
paper, we propose that such cat states can be realized for 



the case of an atomic BEC trapped in a double-well po- 
tential by employing a magnetic Feshbach resonance to 
tune the atomic interaction from repulsive to attractive 
[34H3H . It is found that for suitable parameters, the ini- 
tial state will evolve into a cat-like state after a certain 
period of time, and then it will revive back to the initial 
state. 

Despite these theoretical advances, there remain 
many challenges, some of which we address in the 
present manuscript. One challenge is the problem 
of experimental verification of the coherence between 
macroscopically distinguishable states, which is related 
to the challenge of understanding the role played 
by environmentally-induced decoherence in collapsing 
the cat-state onto a statistical mixture of the two 
macroscopically-distinguishable possibilities. In this pa- 
per we focus in detail on the decoherence which will occur 
if laser-light is used for trapping and/or probing of the 
cat states. In addition the influence of atom-atom in- 
teractions on the proposed detection scheme is studied 
in some detail. As suggested in |36j . we show the re- 
vival of the initial state can be used as an unambiguous 
signal that the supposed cat-state is indeed a coherent 
macroscopic superposition, as opposed to an incoherent 
mixture. We also demonstrate, however, that collisional 
de-phasing can mask the revival without a true loss of 
coherence. 

The effects of decoherence on these cat-like states have 
primarily been discussed for the case of coupling to the 
thermal-cloud of non-condensate atoms. There is dis- 
agreement in the literature regarding whether or not 
this effect is significant [2(| |3(| 13]} • We note, how- 
ever, that these authors treat the non-condensate frac- 
tion as a Markovian reservoir, which implies an infinitely- 
short memory (relative to condensate evolution times) for 
correlations between the condensate and non-condensate 
terms. We question this assumption, primarily due to 
the fact that the non-condensate fraction typically has 
a very low temperature and density. In addition, the 
non-condensate atoms are trapped inside the condensate 
volume by the trapping potential, thus there is no ir- 
reversible loss of information due to propagation of the 
entangled atoms outside of the system volume. Hence, 
we suspect that a non-Markovian treatment is necessary 
to accurately evaluate the thermal-cloud induced deco- 
herence. However, modeling this is a difficult task that 
we do not address at present. Instead, we assume that 
this effect is negligible and concentrate primarily on an 
analysis of a new source of decoherence: that due to 
spontaneous scattering of far-off-resonant photons from 
laser fields which we presume are used in either the trap- 
ping and/or probing of the system. In this case, the 
correlation time of the reservoir is governed by the ir- 
reversible loss of the scattered photons from the system 
volume at the extremely fast rate of L/c where L is the 
atomic mode size and c is the speed of light. Hence, the 
Markov approximation works extremely well for describ- 
ing the physics of photon-scattering. As is well known, 
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this approach leads to effective interactions as well as 
decoherence, effects which include induced dipole-dipole 
interactions and collective spontaneous emission (super- 
radiance) in the atomic ensemble. 

When analyzing the effects of environmentally-induced 
decoherence we are interested in answering two basic 
questions: (f ) What limitations are imposed on the abil- 
ity to create cat-states, and revive the initial condensate 
states by the presence of decoherence, and (2) how long 
can one 'hold' the cat-state by decreasing the tunnel- 
coupling to zero, without collapse of the cat-state onto a 
mixed state. The second question is raised in anticipa- 
tion of applications in precision measurement and quan- 
tum information processing where the cat state is oper- 
ated on and/or held ready for future operations during 
a complicated protocol. We therefore need to determine 
the lifetime of the cat states in the presence of decoher- 
ence, as well as the effect of decoherence on the dynamical 
evolution from the initial state into the cat-state. 

The organization of this article is as follows. In Sec- 
tional the model of bimodal BEC in Double- well system 
is introduced and the cat states in such a system are de- 
fined; in Section lTTTl two schemes of creating cat state are 
discussed: adiabatic evolution of the ground-state and 
and dynamical evolution of a non-stationary state. We 
demonstrate that the former is experimental unfeasible 
due to the fact that the end state is nearly degenerate 
with the first-excited state. We then propose generat- 
ing a cat-like state through dynamic evolution following 
a sudden flipping of the sign of the atomic interaction, 
accomplished via Feshbach resonance. In Section llVl the 
method of detecting the cat-state via revival of the initial 
state is discussed in detail. In Section we investigate 
the decoherence of cat-states and a master equation is 
derived for laser-induced decoherence. In addition, the 
lifetime of cat-states in the presence of laser fields is de- 
termined. In Section IVTI the dephasing effects due to the 
nonlinear two-body interaction is studied. And finally, 
we briefly summarize our results in Section IVIII 



II. THE MODEL 

Our proposed scheme involves loading a BEC into a 
double-well potential and then employing a Feschbach 
resonance to vary the atomic scattering length in order 
to produce a Schrodinger cat-like state. For present we 
consider only the case of a spatial double-well potential. 
However, an analogous system can be formed by using a 
Raman-scheme to couple two hyperfine states in a spinor 
BEC. The primary difference between this model and the 
double- well system lies in the presence of collisional inter- 
actions between the two modes, due to the spatial over- 
lap of the hyperfine modes. Our current model neglects 
these collisions, as they are negligible for the case of well- 
separated potential minima. Most of our results should 
apply to both systems, however, the effects of inter-mode 
collisions will be considered in future work. 



We begin our analysis from the usual many-body 
Hamitonian describing atomic BEC, 

H = J d 3 r^(r)[-^-\7 2 + V trap (r)}^(r) 

+ Y J d 3 r¥ (r)&(r)V(r)$(r), (1) 

where m is the atomic mass, Vt ra p is the trapping poten- 
tial and is the annihilate field operator for atoms in 
Heisenberg picture. The two-body interaction strength 
is given by Uq = Anafi 2 /m, with a denoting the s— wave 
scattering length. In our model, the BEC is assumed to 
occupy a symmetric double-well potential. 

The single-atom ground state and first excited state 
this system we denote as ip s (r) and ipa{r), respectively, 
where tp s (r) is symmetric with respect to the double- well 
symmetry and tp a (r) is antisymmetric. We denote the 
energy gap between these states as hr, which is typically 
much smaller than the gap between the first and second 
excited states. Because of this gap, we assume through- 
out that the atomic population is restricted to these two 
modes alone. It is convenient to introduce two localized 
states Vi, i>R E3> 

Mr) = [<M r ) + V>a(r)] , 

<Mr) = ±= [V> s (r) - ^a(r)] . (2) 

Expanding the field operator ^(r) onto these two modes 
gives 

*(r, t) = c L (t)^ L {T) + & R (t)1> R (T), (3) 

where cl and c R are bosonic atom- annihilation operators 
for the left and right modes, respectively. Inserting this 
expansion back into the many-body Hamiltonian Q , one 
recovers the two-mode version of Bose-Hubbard model 

H = -hT(c[c R + c R c L ) + hg(c^c 2 L + 4 2 c^). (4) 

In deriving this expression, we have implicitly assumed 
that atomic collisions in the overlapping region of the 
two modes are negilibile, valid under the assumption 
that the two modes are well separated. The intra-well 
two-body interaction strength is indicated by g, where 
g = Uq/2% J d 3 r\ip L (r)\ 4 . Note that in the deriving of 
the bimodal Hamiltonian the wavefunction for each 
mode is assumed to be independent of particle number 
in the trap, and terms proportional to the total atom 
number-operator have been dropped, as atom number is 
assumed to be a conserved quantity. 

The quantum state of this bimodal system can be ex- 
pessed as a superposition of Fock-states \n) 

N/2 
n=-JV/2 
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where 



III. GENERATION OF CAT STATES 



(c' L ) N/ 



2+n(J 



)N/2 



y/(N/2 + n)\(N/2-n)\ 



(6) 



Here, N is the total number of atoms in condensate which 
for convenience we take as a even number, and n denotes 
half number difference between the two modes. The two 
corresponding operators are 



N 



c\c L 



\{C ] L C L -$ R C R ). 



(7) 



In this representation, the conventional Schrodinger cat 
state is defined as 



|cat) = -±=(|jV/2) + |-7V/2», 



(8) 



where the N particles have equal probability to be all 
in the left or all in the right mode. This cat-state is an 
ideal maximally- entangled state, however it will be dif- 
ficult to create such a state in this bimodal system. In- 
stead, a more experimental accessible cat-like state is one 
containing two- well-separated wave-packets in the distri- 
bution of c n . For example, a typical cat-like state can be 
of the form 



(« 



c„ ex e 



+ e" 



(9) 



with 7i > a. 

In the present bimodal system the class of cat-like 
states is more general than the Gaussian form @. In- 
stead, a useful cat like state only need to satisfy two 
conditions: (i) the probability distribution in Fock-space 
should be approximately symmetric; and (ii) the two 
wave-packets should be well separated in order to corre- 
spond to macroscopically distinguishable states. In the 
present model, condition (i) is satisfied provided the ini- 
tial state is symmetric, due to the symmetry of Hamilto- 
nian J3J . As we will see is Section |3 decoherence may 
break this symmetry. Nonetheless, to distinguish cat-like 
states from non-cat state, it is convenient to introduce a 
projection operator P cat 



N/2 

j2(\n)(n\ + \-n)(-n\ 



(10) 



This operator acts as a filter which picks up cat-states 
with a minimum wave-packet separation of 2n c . And the 
expectation value (P ca t} can be a efficient way to deter- 
mine whether a state is cat-like. It is straightforward 
to see that for a cat-like state we will have (P ca t) ~ 1, 
whereas for a non-cat state the trace will be noticeably 
less than unity. 



As shown both theoretically and experimentally , the 
s-wave scattering length of cold atoms can be modulated 
by applying a varied magnetic field. In the vicinity of a 
Feshbach resonance, it takes the form 



a = ao(l — 



AB 



), 



(11) 



where ao is the background value, B is the resonant 
value of the magnetic field, and AB is the width of the 
resonance. By tuning the magnetic field B, the scat- 
tering length a can be efficiently tuned in an adequate 
range based on the Feshbach resonance, as demonstrated 
in recent experiments. Thus, the two body interaction 
term g in Hamiltonian , which is determined by scat- 
tering length a, is implicitly assumed to be an arbitrary 
adjustable number. 

Before discussing the methods of generating cat-states, 
we note that attractive BEC's can collapse due to in- 
stability, as been demonstrated in experiments ji^. A 
trapped BEC, however, can be stable against collapse 
under the condition of low atomic density |43j . In terms 
of the BEC atom number N, the attractive interaction 
strength g and the energy gap between the first-excited 
state and ground state Aw, the stability condition can 
be written as 



N < 



Alj 

W 



(12) 



In the present discussion one can increase the energy gap 
between the ground and excited state by adjusting the 
trapping potential, and in this way, we can ensure our 
system will not collapse when the sign of the scattering 
length is switched from positive to negative. 



A. Adiabatic Evolution 

In this subsection, we begin by briefly discussing the 
possibility to create cat-states via adiabatic manipula- 
tion of the many-body ground state [23, , achieved by 
a continuous variation of the interaction strength g from 
a positive to a negative value. In order to better under- 
stand the adiabatic evolution of the ground state, we first 
examine the ground state under certain parameter values 
for which it can be determined exactly. 

The first case we consider is that of zero tunneling, in 
which case the bimodal Hamiltonian can be reduced 
to 



H = 2hgh 



0,2 



(13) 



The ground state of this Hamiltonian is determined solely 
by the parameter g. For g > 0, corresponding to re- 
pulsive interactions, the ground state will be |0) in our 
number-difference representation 10, i.e., an insulator 
state. While in the case of g < 0, the ground state is 



5 



doubly degenerate, corresponding to any arbitrary super- 
position of the two extreme states \N/2) and \—N/2), and 
its orthogonal counterpart. These states can be written 
as 



1 



^a 2 + (3' 



:{a\N/2)+0\-N/2), 



(14) 



with a, f3 being arbitrary complex numbers. Note when 
a = f3, one recovers the Shrodinger cat state ©■ 

In the other extreme case, if there is no atomic in- 
teraction but only the tunneling term, the ground state 
of Hamiltonian JJJ becomes the two-site analogue of the 
superfluid phase |23| 



N/2 

\g) = 2-*V 2 ^ 




In), 



(15) 



which is a binomial distribution of Fock states, and 
thus approximately Gaussian when N is a large num- 
ber. Henceforth wc shall refer to this class of states as a 
'coherent state', meaning that each atom is in a coherent 
superposition of both modes and there are no atom-atom 
correlations or many-body entanglement. 

The more general ground state of our model, with the 
presence of both collisions and tunneling, can be under- 
stood as a smooth cross-over between these two extremes. 
First, under the condition of a repulsive atomic interac- 
tion, the ground states of this system can be well ap- 
proximated by a single Gaussian-like distribution of c n 
peaked at n = 



c n oc e 



7(2° 



(16) 



where the relative atom number dispersion a is obtained 
by minimizing by the mean-energy of Hamiltonian. Sec- 
ond, under the condition of attractive two-body interac- 
tion, the ground state of this system will be split into two 
separated wave-packets. Thus instead of ifTfif . a good ap- 
proximation of this state turns out to be a superposition 
of two Gaussian distributions © 3l], where again the 
minima center hq and spread width a are obtained by 
minimizing the mean-energy. For repulsive atomic inter- 
actions, the energy minima is found at hq = 0, which in- 
dicates that the atom number distribution does not split. 
For attractive interactions, hq ^ 0, and under the condi- 
tion of n > (j, the wavepacket of ground state will split 
and become a superposition of two well-separated com- 
ponents, i.e. it becomes a cat-like state (O- Note for this 
state, the relative phase for any Fock basis pair \n) and 
| — n) is zero. 

Based on these considerations, it is of theoretical in- 
terest to determine whether or not the Shrodinger cat 
state in this bimodal BEC system can be generated by 
an adiabatic process |30L l3l| . As a basic scheme, the sys- 
tem would be first prepared in the many-body ground 
state with positive g, corresponding to a coherent state 
centered at n = 0. Then by slowly tuning g to be nega- 
tive, the distribution will split and finally end up with a 
Shrodinger cat like state under the condition of uq > a. 
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1: The energy gap E gap as a function of ratio g/r. Here, 
gap — Ei — Eq, with Eq, Ei being the eigen energies of the 



ground and first excited state respectively. 



There is a significant difficulty with this proposal, how- 
ever, in that when the two wavepackets in Fock space are 
well separated, the ground state, which is the desired cat 
state, is nearly degenerate with the first excited state, 
as shown in Fig.JTJ. It is observed that when g is pos- 
itive, there is a relatively large energy gap between the 
ground and first excited states. Then as g decreases, the 
energy gap decreases so that at some point the energy 
gap becomes exponentially small. In Fig. [21 we show the 
corresponding wavefunctions of the ground and first ex- 
cited states c(n) with different ratios g/r. It is clear that 
the ground state does not evolve into cat state until the 
bimodal system enters the degenerate regime. In the de- 
generate regime, any small perturbation can mix the two 
states. As the first excited state is nearly identical with 
the ground state, but with a 7r phase-shift between the 
left and right wave-packets, an arbitrary superposition 
of these two states could easily correspond to a localized 
state, which is equivalent to a collapse of the Schrodinger 
cat wavefunction. Strictly speaking, this degeneracy im- 
plies that to adiabatically create the cat state, the process 
requires an exponentially long time. Thus experimental 
realization appears impractical, particularly when con- 
sidering the inevitable exposure of the system to noise. 



B. Dynamic evolution 

An alternate approach to a Schrodinger cat state, 
based on dynamical evolution, has been studied for some 
time [2(| . In this subestion, we show that it is feasible 
to dynamically create the cat state from a coherent state 
by using a Feshbach resonance to make a sudden change 
in the scattering length. Our proposed scheme is as fol- 
lows. First, our system is prepared in the ground state 
with repulsive interactions and strong tunnelling between 
wells. The could be accomplished, e.g. by preparing a 
BEC in a single well, and then adiabatically raising a 




FIG. 2: The wave-packet c(n) of the ground state (Blue solid) 
and first excited sate (Green dashed) with (a) g/r = 0.01, (b) 
g/r = —0.01 and (c) g/r — —0.02 respectively. 



barrier to divide the well into two equal parts. Note 
here, the quantum state of this system is coherent in the 
sense that each atom in the condensate is independently 
in a superposition of the two localized states. The in- 
teraction parameters in this initial stage are chosen to 
satisfy g = g-j^ as this will lead to the optimal cat-like 
state attainable in this scheme |4l|]. Once this state is 
established, the Feshbach resonance is used to achieve a 
sudden switch of the sign of the scattering length, i.e. 
we go from g = to g = — g-^- In order to avoid 
the collapse of the condensate when scattering length a 
becomes negative, the amplitude of atomic interaction g 
should be small. At this point, with tunneling still on, 
the initial state is no longer an eigenstate, thus it will 
start a dynamic evolution under the new Hamilton. 

The evolution of the initial coherent state after the 
change in the sign of g is shown in Figure^ It is seen that 
cat-like states are formed periodically, and between two 
consequent cat states the system approximately revives 
back to the initial coherent state. We have determined 
that the best cat state appears at time t = 14.5/r, which 
corresponds to the final state shown in Fig. [3] Once 
the cat state is formed, one needs to freeze the dynamic 
evolution at a time when the system is in the cat state. 
This can be accomplished by suddenly rising the barrier 
between the two modes to reduce the tunneling coefnent, 
t, to zero. 

Since the present proposal is based on a dynamic pro- 
cedure, predicting the time when the cat state is pro- 
duced t cat is of critical importance. Due to the imposed 
constraint r = 0.9Ng, which is necessary to obtain an 
optimal cat-state, t cat can be expressed solely in terms 



FIG. 3: Evolution of the probability distribution P(n) = |c„| 
versus n, with TV = 100. At t = 0, the initial state is 
taken as the ground state withg = r/(0.97V). Then, the two 
body interaction is suddenly switched to the negative value 
g = — r/(0.9./V). The wavepacket then evolves, where it shows 
a collapse and revival process. It is seen that, at a time of 
t = 14.5/r, a well peaked cat state is obtained with two well 
separated wavepackets centers at around n — 25, 75, respec- 
tively. 



of the tunneling rate r as 



(17) 



where n is determined from our simulations to be k « 
14.5. 

In Figure^] we show the dynamics evolution of projec- 
tion operator P ca t- Again, it is shown that the optimum 
cat-like state forms at around t = t ca t = 14. It and reap- 
pears again at intervals of t = 10/r. It is also noticed 
that between each cat state, there are valleys of P ca t, 
which correspond to revivals of the initial coherent state. 

It is also of interest to study the relative phases be- 
tween the Fock states in the cat-like state generated in 
this manner. The symmetry between left and right modes 
is automatically guaranteed in our model since our sys- 
tem Hamiltonian (0J and our initial state are both sym- 
metric. In Figure [SJ we see that the phase of created cat 
state is exactly symmetric around n = 0. In addition, we 
see that the phase flattens somewhat in the vicinity of 
n = ±25, where the peaks in the probability are located. 

We have claimed that imposing the condition r = 
0.9Ng is necessary to produce the optimal cat state. It 
is worth pointing out that with other choices, left-right 
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FIG. 4: Time evolution of the expectation value of the pro- 
jector Peat with system parameters chosen as in Fig. © and 
a filtering window of width n c = 15. A well peaked cat state 

is created at t = 14.5/r, 24/r, 34/r Note the P cat reaches 

a minimum at t » 10/r, 20/r, 30/r..., which corresponds to 
the revival of the initial coherent state. 
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FIG. 6: The measurement P ca t versus interaction strength 
variation, with g — Q ^ N (1 + 5 B ). Other system parameters 
and projection operator is identical to Fig.l^IJ. It is shown 
that in the vicinity of 5 g = 0, the cat state is created. Note 
that the projection is measured at different times for each S g , 
so that the value of P ca t is maximized. 
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FIG. 5: Phase distribution of the created Schrodinger cat 
state in Fock space. The y axis is the phase 6 n = Arg{c n ). 



symmetric states can be created, but these are neither 
well-separated nor double peaked. This is shown in Fig- 
urelHl where we plot the expectation value of the projec- 
tor P cat versus S g , where we have taken g — (1 + S g ). 
In this figure, we see the cat-like state is only created 
in the vicinity of S g = with an uncertainty of approx- 
imately ±5%. That means to create a cat state in the 
present system, a precise control of the atomic interaction 
g or tunneling rate r is required. It is important to note 
that in generating Figure [S] we have varied the evolution 
time t cat for each parameter choice, in order to maximize 
the projection P ca t- Thus Figure can be interpreted as 
verifying the conditions g — and t cat = 14.7/r. 

Up to now we have only considered the ideal situation 
where the initial state is exactly prepared in the ground 
state and the atomic scattering length can be tuned in 



a precise way. This may not necessarily hold true un- 
der experimental conditions, where random fluctuations 
in system parameters may affect the formation of the 
cat-state. Here, to test the susceptability of the present 
proposal to such fluctuations, we study the dependence 
of final cat-state on the precision of the initial state by 
adding random noise to the exact ground state. Thus, 
the initial state becomes c n + 5 r . The random noise, 
5 r , is a complex with an amplitude drawn from a Gaus- 
sian distribution, and a completely random phase. The 
probability distributions for the amplitude and phase are 
therefore 

P(\S r \) = -i-e-MV-?, 
V 7ro > 

P(0) = (18) 

Figure shows the resulting cat-state as a function of 
oy . We see that a well peaked cat state can be produced 
up to ay « 0.05, which is about co(t = 0)/6, i.e, the 
peak of initial coherent state we make the cat state from. 
Thus this method of generating a cat state appears robust 
against noise in the initial state, which could be do to a 
finite temperature. 

Another possible source of error in creating the cat 
state is from the inaccuracy in control of the atomic in- 
teraction strength g and/or the tunneling rate r. This 
effect is simulated by adding a Gaussian-distributed ran- 
dom number to negative g. We note that this is not 
equivalent to time-dependent fluctuations in g and r, but 
rather represents imprecision in the control over the val- 
ues of the parameters. The result is shown in Fig. |5J 
where it is seen that the evolution system is sensitive 
to the noise of the atomic interaction g, with only about 
±2% deviation being tolerable. However, it does not nec- 
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FIG. 7: The probability distribution P(n) at time t — t ca t 
is plotted versus the level of random noise in the initial state 
ay. We see that even at ay = 0.05 the output state is still a 
good cat-state. This corresponds to random fluctuations in 
the initial c n 's at about 5% of their exact ground state values 
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n 



FIG. 8: The phase distribution 8 n — arg(e n ) at t = t cat is 
plotted for the case of random noise in the initial state with 
o r = 0.05. We see that the phase distribution agrees well 
with Fig. |S] up to a non-physical overall phase-shift. The 
disagreement occurs only in the regions where the amplitude 
of Cn is small. 



essarily mean that one can not create a cat-state without 
the well defined g. Instead, when g is shifted, the time 
need to create cat state t cat is also shifted, and one can 
still obtain the cat-state on another time. The primary 
difficulty here is that if the variation in g is unknown for 
a given run, then the precise time to stop the dynamics, 
t ca t is not known a-priori. 



IV. DETECTION OF CAT STATE 

Presuming a cat-like state can be created in the above 
manner, it remains a challenge to experimentally ver- 




FIG. 9: The effect of imprecise control over the atomic inter- 
action strencth g. The noise of g is introduced by adding a real 
random number, where a r is the width of the corresponding 
Gaussian distribution, to the value of g in Hamiltonian ||1J. 
The cat state is measured at a time of t — 14.5/r and plotted 
as a function of a r . It is shown that uncertanties in g up to 
3% can be tolerated. 



ify the coherence between the macroscopically distin- 
guishable states. Neverthless, a possible method to de- 
tect Schrodinger cat states in this bimodal BEC's lies 
in detecting the revival of the initial coherent state. In 
Fig. ^1] we show the evolution of the cat state starting 
from t = t ca t with the system parameters still satisfying 
g = q N . This shows that under the condition of no 
decoherence, the system will evolve back to a coherent 
state after a certain period of time, which we denote as 
t rev . The revival time t rev again is estimated from our 
simulations, yielding t rev = 7/r. 

In the presence of decoherence, however, the system 
may not be exhibit a distinct revival of the coherent state. 
To illustrate this we first consider two states which we 
define as the partially incoherent state and the totally 
incoherent state. These states are defined by the intro- 
duction of a coherence 'length' in Fock-space, I co h, so 
that coherent superpositions between Fock-states \n) and 
|m) are destroyed for \n — m\ ^ I C oh, while superposi- 
tions satisfying \n — m\ < l C oh are relatively unaffected 
by decoherence. The partially incoherent state is then 
defined in the sense that the coherence length I co h is 
small compared to the distance between the two sepa- 
rated peaks of the cat state, but larger than the width of 
either wavepacket, i.e., I co h < no yet I co h > a in Eq. 10. 
For the completely incoherent state, on the other hand, 
the cat state is taken to have collapsed to a single Fock 
state, i.e. a mixed state maintaining the probability dis- 
tribution of the cat-state but with to coherence between 
Fock-states. 

The evolution of those two states in time is shown in 
Fig. Illland Fig ll2l respectively. Unlike the cat state, the 
incoherent states will not exhibit a revival of the initial 
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FIG. 10: Revival dynamics from the cat state back to the 
coherent state. This figure shows the continued dynamical 
evolution following that shown in Fig[3]with g — —0.01. At 
t — 14.5/t, the quantum system is in a cat state, then at a 
time of t » 21, the initial coherent state is revived. 



coherent state at t = t rev . The partially incoherent state 
evolves instead into a disinct three-peaked state, whereas 
the completely incoherent state exhibits small oscillations 
with no discernable peaks after a very short collapse time. 

Based on these considerations, one prospective method 
to verify the coherence of the cat state would involve: i) 
first, measuring the probability distribution, P(n), which 
is accomplished by repeated measurements of the atom- 
number in each mode; ii) then after forming the cat-state, 
the state should be held at r = for a time t = thoid an d 
then with tunneling restored, the revival of the initial co- 
herent state should be observed after a time t rev ; iii) a 
known source of decoherence should be added while hold- 
ing the system at r = for duration t^oWi after which 
a lack of revival should be observed. To add decoher- 
ence, one could simply add a laser field to the system. 
This will introduce decoherence due to spontaneous pho- 
ton scattering, as described in Section If the revival 
of the initial state can be observed in the case of no ad- 
ditional decoherence, but disappears when the decoher- 
ence is added in a controlled way, than one should have 
a strong claim for the observation of coherence between 
macrscopically distinguishable states. 

V. LASER-INDUCED DECOHERENCE AND 
LOSSES 

While trapped ultracold atoms are typically well- 
isolated from the environment, the presence of lasers 
for trapping and/or detection of atoms will polarize the 
atoms and thus couple them to the vacuum modes of 
the electromagnetic field. This coupling allows for spon- 
taneous Rayleigh scattering and/or diffraction of laser 
photons, which can lead to losses and/or decoherence of 




FIG. 11: Dynamical evolution of a partially incoherent cat 
state. The parameters are chosen the same as Fig llOl At 
t — 14.5/t, the quantum system is in a mixed state with no 
coherence between the two wavepackets. Then as it evolves, 
the system exhibits oscillations, but will not evolve back to 
the initial state. 




FIG. 12: Dynamic evolution of completely decoherent cat 
state. At t = 14.5/t, we let the system collapse onto a sta- 
tistical mixture of Fock-states. It is seen that the dynamic is 
completely disordered. 



the atomic system. If the scattered light field acquires 
sufficient information to determine the atom-number dif- 
ference between the left and right wells, then this spon- 
taneous scattering will lead to a dynamical collapse from 
the cat-state into a mixture of Fock states, an example 
of environmentally-induced decoherence. 

In addition, the condensates are in contact with a ther- 
mal cloud on non-condensed atoms. In |36| . it was deter- 
mined that such a coupling will lead to a rapid collapse of 
the cat-state. We note, however, that the treatment in 
|36j | made the approximation that correlations between 
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the condensate and the cloud decay much faster than the 
time-scale of the condensate dynamics, i.e. the thermal 
cloud was treated as a Markovian reservoir. At present 
we question the validity of the Markov approximation for 
such a system. The coherence time can be quite long be- 
cause (i) the cloud temperature is very low and the cloud 
is very dilute (so that collisions between cloud atoms are 
rare), and (ii) the cloud is confined, so that information 
does not propagate away from the condensate region. A 
non-Markovian treatment of thermal-cloud-induced de- 
coherence is beyond the scope of this work. At present, 
we focus primarily on the decoherence induced by cou- 
pling to the electromagnetic field. 

The structure of this section is as follows: in subsec- 
tion A, we will first investigate the decoherence mechan- 
ics of condensates induced by laser scattering, where the 
general master equation is derived. Then, in subsection 
B, we discuss the condensate losses due to the inelastic 
collision, and the loss rate is obtained for one-mode con- 
densate. In subsection C, we show the decoherence of 
the two-mode cat states in the elastic scattering regime, 
defined as the regime where the photon-recoil is not suffi- 
ciently strong to remove the atom from its initial center- 
of-mass state. The system dynamics under this loss-less 
decoherence is studied in subsection D.. Lastly, the atom- 
loss and decoherence of cat-states in the inelastic regime 
is studied in subsection E, where we find that the cat- 
states created via dynamic evolution are less susceptible 
to decoherence-induced collapse than the extreme cat- 
state, i.e. the true ground-state of the system, which we 
have argued is inaccessible via adiabatic evolution. 



A. Master equation 

If there is an optical component to the trapping poten- 
tial, the atoms in the BEC will have an induced electric 
dipole-moment which oscillates at the laser frequency. 
This dipole-moment will couple the atoms to the elec- 
tromagnetic vacuum, resulting in spontaneous Rayleigh 
scattering and/or diffraction of the laser beam. In ad- 
dition, there will be laser-induced many-atom effects, 
primarily dipole-dipole interactions and collective effects 
such as superradiance. As we will see, these effects will 
lead to decoherence and collapse of the cat-states, even 
in the absence of recoil-induced losses from the two cat 
modes. This decoherence is associated with the possibil- 
ity for scattered and/or diffracted photons to carry away 
information sufficient to reveal the distribution of atoms 
between the two modes. 

In order to model these effects we need to derive the 
master-equation which governs the density operator of 
the atomic field. In the presence of laser coupling to a 
BEC system, the total Hamiltonian for the system, which 
includes the Hamiltonian for the atoms in condensate H s , 
the reservoir of vacuum electromagnetic field H s and the 
interaction between them in the dipole approximation 



V sr can be written as 

H = H s + H r + V sr , 



(19) 



H s = J d 3 rA*t(r)* e (r) + J d 3 rft(r)^(r)* g (r), (20) 



E 

feA 



(21) 



V„ 



E 

feA 



(22) 

where x & e (r) and ^(r) are the annihilation operators for 
atoms in the excited and ground state respectively, and 
SkA is the photon annihilation operator for momentum 
Kk and polarization state A. In addition, uj ei lul and LUk 
denote the frequency of atomic excited state, laser field, 
and the reservoir photons respectively, while A = w e -WL 
is the detunning of laser field from the excited state. The 
atom-laser interaction is governed by the Rabi-frequcncy 
f2(r), while the atom-reservoir interaction is governed by 

9k\ -- 



2 ^v~& ■ e/cA, with d being the atomic dipole 

moment in mks units. We note that for the purpose of 
deriving the master equation it is not necessary to include 
terms in H s governing the free evolution of the ground 
state, as this involves very slow time scales relative to 
the dynamics of the excited-state/EM- vacuum system. 
These terms can be added to the system Hamiltonian 
at the end of the calculation, as any energy-shifts they 
introduce will be negligible compared to the natural line- 
width of the excited state. 

The quantum state of the system+reservoir is de- 
scribed by the density operator p sr . An effective equa- 
tion of motion for the reduced system density operator, 
p s = tr r {p sr }, can be derived via second-order perturba- 
tion theory by following the standard approach (see for 
example [4J|). For the case of a zero-temperature bath 
this leads directly to the master equation 



p s (t) = -iH s p s (t) 



d 3 rd 3 r' 



L(r - r')¥ e (v)* g (rM(r')4? e (r')p s (t) 



L*(r - r')*U r )*e(r)p s (0*l(r')* s (r') 



H.c. 



(23) 



where 



L(r-r') = Y^\9k\\-' 

feA 



2 tk-(r-r') 



{ir5(uj k -ul) + ■ 



iV 



with V indicating a principal value. In the limit of in- 
finite quantization volume, the summation becomes an 
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integral, which can be done analytically |4^, yielding 



L(r-r') 



3r 



C 



(1 -3 cos 2 



1 



c 2 



where T = lo\<P / (Sireofic 3 ) is the single-atom sponta- 
neous emission rate, C = ki,\r — r'\ and is the angle 
between the dipole moment d and r — r'. The real part 
of L(r — r') will contribute to decoherence in the master 
equation I23|) . while the imaginary part contributes an 
energy shift due to photon exchange between atoms. In 
the limit |r — r'l — > we find 



L(0) = -+iS, 



(26) 



where the imaginary part, S is an infinite quantity, re- 
flecting the fact that the standard two-level atom-field 
interaction model is not properly renomalized. Physi- 
cally, this term is a Lamb-type shift in the energy of the 
excited-state due to interaction with the vacuum modes 
of the EM field. In our system we can always choose a 
rotating frame to absorb this shift, so that we can take 
(5 = without loss of generality. 

For most optical traps, the laser is detuned very far 
from the atomic resonance frequency, so that the excited 
state occupation number is <C 1. In this regime it is 
useful to perform an adiabatic elimination of the excited- 
state operators. The resulting master equation, involving 
only ground-state operators, is 



p s {t) = -iH s p s (t)-i / d 3 r 



|0(r) 



-\&t(r)*(i 



A 



d 3 rd 3 r / n(r) L<yr A9 r ^ n*(r')^(r)^{r)¥(r')^(r')p s {t) 



n(i^)tft(r)tf (r)p,(t)¥t(r')tf (0 + H.c, 
I 



(27) 



where ^(r) = ^ g (r). The details of this derivation are 
presented in Appendix IVIIII 



B. Condensate losses: elastic and inelastic regimes 

We will now employ this master equation to describe 
decoherence and/or atom-loss in our two- mode atomic 
system. Underlying the dynamics governed by master 
equation 127(1 is the exchange of photons between atoms 
and between an atom and the EM vacuum. As such 
photon exchanges can involve significant momentum ex- 
changes due to photon recoil, it is useful to first de- 
fine elastic and inelastic regimes with respect to the ini- 
tial atomic center-of-mass state. If the spatial size of 
the initial atomic mode is small compared to the laser 
wavelength, then one would be in the elastic regime, as 
photon-recoil would not be sufficient to carry the atom 
outside of the momentum distribution of the initial state. 
In the inelastic regime, on the other hand, the initial 
mode is large compared to the laser wavelength, which 
implies that it is narrow in momentum space compared 
to the photon momentum, so that the recoil from a sin- 
gle photon is sufficient to remove an atom from its initial 
mode and place it in an orthogonal mode. 

In the strongly elastic regime, there will be negligible 
loss of atoms from their initial modes. As we shall demon- 
strate, the master equation 1(27(1 may still lead to decoher- 
ence/collapse of the two- mode cat-states. In the inelastic 
regime, on the other hand, the predominant effect of the 
atom-field interaction will be scattering of atoms out of 
the initial modes and into a quasi-continuum of modes. 
This loss will be accompanied by a collapse of state of 



the remaining atoms into a mixture of states, which may 
or may not be cat-like in themselves. The details of this 
decoherence process are given in subsection IV El 

To derive the atom decay rate of a single atomic-field 
mode, we will need to make use of the commutation re- 
lation 



[c t c,* t (r)*(r)] 



c)c^{r)-<j)*(r)¥(r)c, (28) 



where c is the annihilation operator for atoms in mode 
</>(r). Then, from the master equation ((27(1 . it is straight 
forward to obtain the evolution of the mean mode occu- 
pation no 



= (etc) as 



n 



<P*(r)5 3 (r-r')(c^(r')) 



(29) 



where G(r,r') = fi(r)£(r — r')fT(r'). Making a single 
mode approximation, 



#(r) 



r c, 



(30) 



one obtains 

n = - 



2 

A2 



d 3 rd 3 r' 31 [G(r, r')] |^(r)| 2 

x [5 3 (r-r')-|^(r')| 2 ] no. (31) 



To better understand the physics described by this 
equation, we make the simplifying assumption f2(r) ss 
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FIG. 13: Dependence of the inelastic loss parameter £ on the 
spatial size of condensate, which is indicated by the dimen- 
sionless parameter a as shown in Eq. 13411 . It is seen that £ is 
approximately unity when the size of condensate is much less 
than the laser wavelength, i.e. a <C 1, which means there is 
negligible atom loss. As the condensate size approaches half 
Al, £ decreases to around 0, corresponding to the maximum 
loss rate. 



ing k L ■ d 



where k^ is the laser wave- vector, satisfy- 
0. This leads to 



n 



^r(i 

A 2 1 



n o, 



(32) 



where 



d 3 rd 3 r'n[e 



ik.L-(r— r') 



L(r-r')]|0 fl (r)| 2 |0 g (r')| 5 



(33) 

In the elastic regime we have |</>(r)| 2 ~ S(r) relative to 
L(r — r'), which leads to £ = 1 so that there is no atomic 
loss, as expected. As the size of condensate grows, £ 
will decrease, as shown in Fig. 1131 which plots £ versus 
condensate size for the case of a spherically symmetric 
Gaussian condensate wavefuction 



1 



(CTA L ) 3 / 2 7rV 4 



(34) 



where a is the ratio of the condensate size to the laser 
wavelength, A^. From the figure, we see that the de- 
cay rate of the condensate can be neglected under the 
condition it C 1, while for a mode whose size is com- 
parable to or larger than the laser wavelength, we have 
£ — * so that the population will decay at the expected 
rate |f2| 2 r/A 2 , which is simply the electronically-excited 
state fraction times the excited state lifetime. 



tween the wells. Coupling of this system to the EM vac- 
uum modes will result in decoherence only when scat- 
tered photons carry away sufficient information about the 
number of atoms in each well. It thus follows that if the 
separation between the two condensates, s, is less than 
the laser wavelength A^, no decoherence should occur, as 
the information contained in the photons is diffraction- 
limited to a resolution of no better than a wavelength. 
On the other hand, if s is larger than Al , the locations of 
the atoms can in-principle be determined from the phase- 
information contained in the scattered photons, this ir- 
reversible transfer of information to the environment will 
be reflected in decoherence of fock-space-superpositions 
of different atom-number distributions. 

In the deeply elastic regime, the size of the left and 
right modes is the smallest length scale in the problem, 
so that we can replace the atomic density distributions 
with delta-functions, or equivalently we take 



*(r) = a 5(r 



1)cl+ JS(r+^)c R , (35) 



where cl and c R are the annihilation operators for the 
left and right mode respectively. Inserting this expansion 
into the master equation l|27|l and assuming again for 
simplicity fi(r) = f2e 4ki ' r , gives 



Ps 



2A 2 



& h C L & L C L ps 



CrCRCrCRPs 



+2cos(k L • s)p(s)c[c L c R c R p s - c> L c L p s c f L c L 
-c R c R p s c R c R - cos(k L • s)p,(s)c[c L p 3 c R c R 
- cos(k L • s)p(s)c R c R p s c\c L ) + H.c (36) 

where again the imaginary parts are absorbed via a trans- 
formation to the appropriate rotating frame. The result- 
ing decoherence is governed by the parameter p(s), de- 
fined as 



p(s) ee 2R[L(s)]/r 



= - sin 

2 L 



sin( , 2 wCosC sinC, 

,— h (1 — 3cos S )(- 



c 2 



(37) 



where £ = 2-ks/Xl and 6 S is the angle between the in- 
duced polarization direction d and the relative coordinate 
s. 

From this master equation H3bf) . we derive the equa- 
tion of motion for the fock-space matrix elements p nm — 
{n\p s \m), yielding 



Pn 



IO|2 

r LJ- ( n - m) 2 (l - cos(k L • s)p(s)) Pnm . (38) 



C. Decoherence of two-mode cat states in elastic 
regime 

In the present scheme, N atoms are trapped in a 
double-well potential with a finite spatial separation be- 



This shows that the diagonal matrix elements (m=n) do 
not decay, reflecting the absence of condensate losses 
in the elastic regime. The off-diagonal elements, on 
the other hand, decay as at a rate proportional to 
(1 — cos(k£ ■ s)/x(s))(n — m) 2 . 
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The dependence of fi(s) on the mode separation is 
shown in Fig. O which plots /i(s) versus s/Xl and 
6 S = args. We see that under the condition of s <C Xl 
we have cos(k • s)/i(s) w 1, and therefore no decoher- 
ence. As s increases, /i(s) will decrease and the deco- 
herence rate will increase. In the limit of s 3> Xl, we 
have /x(s) « 0, and the decoherence converges to the rate 

Inm — -%rF(n — m) 2 . It is interesting to note that /x(s) 
is angle-dependent, as seen in FiglTH 

To understand the reason for this decoherence rate (|38|) 
we will consider the limiting cases s 3> Xl and s < Al, 
corresponding to a mode separation much larger than or 
much smaller than the laser wavelength, respectively. For 
the case where the mode separation is large compared 
to the laser wavelength, we have /Lt(s) — > so that the 

decoherence rate is j nm — ^grT(n — m) 2 . Our goal is 
thus to explain the peculiar (n — m) 2 dependence in this 
expression. We begin by noting that in this regime the 
phase information in the scattered photon wavefronts can 
in principle determine whether the photon scattered from 
the L or R mode. Thus the question becomes: on what 
timescale, T, does the radiation held scattered by L or R 
acquire sufficient information to distinguish the state |n) 
from the state |m). This timescale will then determine 
the decay rate for the off-diagonal elements p nm , which 
give the degree of coherence between the states \n) and 
m) . The rate at which photons will scatter from the left 
In 1 2 

mode is ^"Tn 2 , where ul is the number of atoms in the 
left mode. The reason for the n\ dependence, rather than 
the usual ul, is due to Bose stimulation, as the atoms 
remain in their initial mode after elastic scattering. 

The actual scattering of photons is a random pro- 
cess, so that the number of L scattered photons n p 
over a time interval t is drawn from a distribution 
function PL(n p ,riL,t). It is reasonable to assume that 
Pl(ti p , Ul, t) is at least approximately Poissonian, with a 

center at n p (riL, i) = ^p-Tn 2 i, and a width An p (riL, t) = 

^/ n p( n L,t) = ^VTtriL. Based on this assumption, the 
criterion for distinguishing state \n) from state |m) is that 
the distributions pl (n p , til, t) and pL(n p ,mL,t) should 
be distinguishable, so that the scattered photon number 
can be attributed to one distribution or the other. In this 
expression til is the number of atoms in mode L for state 
\n) and This requires minimal overlap between the two 
distributions, which can be approximately formulated as 

n p {n Ll T) + An p (n L ,T) = n p (m L ,T) - An p (m L ,T), 

(39) 

where we have temporarily assumed tul > til- Solving 
this equation for the decoherence rate 7„ m = 1/td yields 



A 2 [niL+nL] 2 
From eq. J7J) it follows that 



(40) 



N 

n L = T -n 



N 

n R = — 



n. 



(41) 



which also hold for m l and m,R with a simple substitution 
n — > m. Inserting this into l|4(J|) yields 



M 2 

A 2 



(42) 



We note that the result for mode R instead of L can be 
found via n <-> —n and m <-> — m, while the result for 
mL < hl can be found via n <-> m, all of which leave the 
result (|42() unchanged. Thus we have verified our inter- 
pretation that the decay of the p n m coherence is governed 
by the timescale on which the information carried by the 
scattered light becomes sufficient to distinguish between 
the states \n) and \m). 

If we consider the opposite case s <C Xl we can no 
longer assume that scattered photons carry information 
concerning which mode they scattered from, due to the 
standard diffraction limit. In this case one might be 
tempted to assume that the total scattering rate is pro- 
portional to n\ + n\ — N 2 /2 + 2n 2 , as photons must 
still scatter from one mode or the other and each mode 
will experience Bose-stimulation of the recoiling atom. 
This would imply that the total scattering rate would be 
different for the states \n) and |m), so that given suf- 
ficient time the environment could learn which state is 
scattering the light and thus destroy the coherence. This 
effect, however, is compensated for by superradiance 0] 
so that the scattering rate is proportional simply to iV 2 , 
which is the same for both |n) and |m). Superradiance is 
a many-body quantum-interference effect whereby atoms 
within a distance less than the optical wavelength experi- 
ence enhanced decay, identical to Bose stimulation, even 
for distinguishable atoms. 

To verify this interpretation one can calculate the scat- 
tered light intensity and see precisely how it scales with 
Ul and tir. This can be accomplished by using Eqs H21|) 
and H22f> to obtain the Heisenberg equation of motion for 



dt 



_ ^ 5 * A e- 4 ^* / d 3 re i ( kl - k > r l' t (r)l(r|,43) 



where we have substituted the adiabatic solution tp e {^) = 
■^^(r) • Substituting the mode expansion (|35fl and not- 
ing that fiL and fin are constants of motion in the absence 
of tunneling gives upon time integration 



fikA(i) = a k A(0)e-^* - t ^ x e-^ + ^ 2 



sin[(u fc - m L )t/2] 
(LOk - u>l)/2 



ei (k-k i ). f ~ i + e - l( k-k L ).|- fl 



(44) 



From this expression we can compute the mean number 
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FIG. 14: The decoherence parameter fi(s) as a function of 8 S 
and s/Xl- From the diagram, /i(s) is approximately 1 when 
the separation of two modes s is much smaller than the laser 
wavelength Xl- As s increases, fi(s) decreases at a speed of 
1/s and eventually approaches zero as s goes much larger than 
Xl- Note that fj,(s) is also dependent on 9 a , where at large s, 
/j,(s) ~ sin 2 8 B . 



of scattered photons via n p (t) = X^kA^icA^k-O' yielding 

n p( f ) = ^2" r + u2r + 2cos ( k i ' s )m( s ) ninfl] i 

. (45) 

where we have converted the sum to an integral in the 
limit of infinite quantization volume and made the ap- 
proximation sin 2 (xt)/x 2 « TrtS(x) when integrating over 
k. This expression reveals an additional component to 
the scattering rate proportional to cos(k£ • s)/i(s)n_L7iij. 
This term is only significant when /it(s) ~ 1, which re- 
quires that the mode separation be comparable to or less 
than the optical wavelength. In addition it depends on 
the exact spacing and orientation of the modes with re- 
spect to the polarization of the lasers, via the implied de- 
pendence of /i(s) on s. We interpret this additional scat- 
tering as superradiance, as it arises from the imaginary 
part of the two-body dipole-dipole interaction. We see 
that in the limit as s/Xl — > we have cos(ki >s)/x(s) — > 1 
so that the total scattering rate is proportional to TV 2 . 
As stated previously, this means that the radiation field 
never acquires sufficient information to distinguish the 
states |n) and |m), thereforet no decoherence occurs. 

D. Dynamics under decoherence in elastic regime 

In the previous section, we derived the master equation 
of the two-mode system in the presence of laser-induced 
atomic polarization. In this section we use numerical sim- 
ulations to observe how the decoherence will affect the 
tunneling dynamics of our system and determine what 
level of decoherence can be tolerated during the dynami- 



cal generation of a cat state. In this section we will only 
focus on the decoherent dynamics in this elastic regime 
where the total atom number can be treated as a con- 
stant of motion. In addition, we will assume that the 
separation of condensates are much larger than the laser 
wavelength, so that the dynamics is determined Eq. i|38|) 
but with /i s = 0. Under those assumptions, it is straight- 
forward to guess the coherence time t co h as 

tcoh ~ iv 2 |0| 2 r' (46) 

by observing that thepeaks of the dynamically created 
cat-state are separated by no = ±AT/4. This is the time- 
scale on which coherence between states separated by 
Ay 4 should decay based on the decoherence rate t4*2)l . 
During the generation of the cat state, it is expected that 
if the coherence time t co h is large compared to t cat , the 
cat state can still be produced. Otherwise, the system 
should evolve into a mixed-state rather than a pure cat- 
state. 

The effects of decoherence on the dynamical evolution 
of cat-states is shown in Figures 1151171 In Fig. ^| we 
show the magnitude of the density-matrix elements for 
the initial coherent state at t = 0. Under the condition 
of no decoherence, this state evolves into a cat-state at 
t = t cat , which is shown in Fig. In Figures IT7I (a) 
and (b) we show the resulting state for the cases of weak 
and strong decoherence, respectively. In Fig. I17f a) the 
decoherence strength is chosen such that t co h — 5t ca t, 
while for the strong decoherence of Fig. 1171 (b). we have 
t C oh = t ca t- H is seen that with t co h = 5t ca t we can still 
create cat-like states, while as the decoherence strength 
increases up to t co h = t ca t, the outcome state is a mixtcd 
state with no coherence between macroscopically distin- 
guishable distributiions. Therefore, to successfully pro- 
duce cat state, the coherence time t co h should be much 
longer than the time need to produce cat state t ca t- 

It is also quite interesting to observe from Fig. 1 171 that 
not only the off diagonal, but also the diagonal of den- 
sity matrix of the final states are affected by decoherence 
effect. In Figure ITS1 we plot the probability distribution 
P(n) = Pnn of several resulting states corresponding to 
different coherence times. In the figure we see that the 
probability distribution of the mixed states formed under 
strong decoherence iis significantly different from that of 
the pure cat state. In other words, if decoherence dom- 
inates the system dynamics, it will not yield the same 
P(n). Thus the probability distribution alone supplies 
a direct way to verify the existence of cat state at time 
tcat, without the need to detect revivals. We note that 
we have so far only established this result for the elas- 
tic regime. We suspect that the elastic regime will be 
much more difficult to achieve experimentally than the 
inelastic regime, due to the necessity of sub-wavelength 
confinement. 

In Figure El we examine the effects of decoherence on 
the revival process, under the assumption that t ca t < 
t C oh, but assuming that the system is held at r = for 
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FIG. 15: Magnitude of the density matrix elements for the 
initial coherent state. 
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FIG. 16: Magnitude of the density matrix elements of the 
dynamical created cat state at t — t cat with no decoherence. 



a time thoidi which may be longer than the coherence 
time. These figures show the magnitude of the density- 
matrix elements at time t — t cat + t} lo id + t rev , calculated 
by direct numerical solution of the master equation (|36[) . 
In Figure EH ( a ) we show the magnitude of the density- 
matrix elements for the case t^oid — and for t co h = oo, 
i.e. no decoherence is present. Note that this figure cor- 
responds to t = 21 in Fig. In Figure H5I (b) we plot 
the magnitude of the resulting density matrix for the case 
thoid — t co h — lOtcot' This choice gives sufficient time to 
destroy coherence between fock states where n differs by 
N/2, corresponding to the distance between peaks in the 
initial cat state. This shows that the revival is effectively 
destroyed by decoherence on the the timescale governed 
by the coherence time Q46p. This figure corresponds to 
the 'partially incoherent state' shown in Fig. 1111 at time 
t = 21, in the sense that the choice thoid — tcoh im- 
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(a) 



lp I 




(b) 

FIG. 17: Magnitude of density matrix elements of the result- 
ing state under weak decoherence (a) and strong decoherence 
(b). System parameters are chosen the same as in Fig.Bl. 
and the resulting states (a) (b) are measured at a time of 
t cat = 14.5/r. 

plies the condition a < I co h < no- Lastly, Figure H5I (c) 
shows the case t hold = \t co h with t co h = Wt cat which 
gives sufficient time to destroy coherence between any 
fock states. We see this case that the system has col- 
lapsed onto an almost pure mixture of fock-states, and it 
is clear that no revival occurs. This figure corresponds 
to the totally incoherent state of Fig. ITS! 



E. Loss and Decoherence in the inelastic regime 

In this subsection we turn to the inelastic regime, 
which takes into account the finite size of the atomic 
modes, so that the scattering of laser photons is not en- 
tirely elastic due to photon recoil effects. The resulting 
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FIG. 18: Probability distribution P(n) of the resulting 
states corresponding to different coherence times. Pa- 
rameters are choose as in Fig. I17L while from 
line (1) to (7), the coherent time t co h is chosen as 
oo, St cat , 3t cat /2, t C at, 3t ca t /4, 3t cat /5, t ca t/2, respectively. 



losses from the two initially occupied modes should nec- 
essarily be accompanied by decoherence, which may or 
may not be sufficient to 'collapse' the two-mode cat state. 
We focus solely on the case where the mode separation is 
large compared to the wavelength, so that we can safely 
ignore the cross terms in Eq. l|36fl . The master equation 
in this situation is found by substituting the two-mode 
expansion ^(r) = 0i(r)cL + 4>b.{^)cr, yielding 



Ps(t) 



[Np s (t)+Z[n 2 L +n R -N]p s (t) 

He, (47) 



2A 2 

£,[n L p s (t)n L + n R p s (t)n R }] 



with £ defined in Eq. (J2HJ). Deriving the equation of 
motion for the density matrix elements is then straight- 
forward, giving 



Pnm (^) 



-T^[(l-ON + ^n-m) 2 ]p nm (t). (48) 



Under the condition of spatially broad (relative to Al) 
condensates, £ ss 0, as seen if Figure ll3l and this equation 
is reduced to 



Pnm (^) 



-r-^2-ivp nm (t), 



(49) 



which shows that the decoherence rate is exactly equal to 
the atomic loss rate (|32fl . The reason for this is that the 
loss of an atom takes us from the manifold with N atoms 
to that with N —1 atoms. This is equivalent torn n± 1 
where the choice of ± depends on whether the atom was 
lost from the L or R mode (- for L and + for R) . Thus one 
might suspect that a careful treatment of the scattering 
of the atom into the quasi-continuum of recoil modes 
would reveal that the loss of coherence on the matrix 
element p nm could be accompanied by a corresponding 
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(c) 

FIG. 19: Effects of decoherence on the revival of the initial 
coherent state. Shown is the magnitude of the density matrix 
elements of the system at time t = t ca t + thoid + trev under no 
decoherence (a), the decoherent resulting state with t^oid = 
tcoh (b), and thoid = N 2 t co h/^ (c). System parameters are 
chosen the same as in Fig.©. 
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appearance of coherence on the matrix element p n ±i,m±i 
which lies in the N — 1 manifold. 

A complete quantum description of this inelastic pro- 
cess should include both the initial cat-state modes as 
well as a quasi- continuum of recoil modes. We can there- 
fore describe the initial unscattered systems with a den- 
sity of product p^ r = p N (g>p r , where p™ is the full density 
operator having m atoms in the initial modes and N — m 
atoms in the quasi-continuum, and p m and p^~ m are 
the corresponding reduced density matrices for the sys- 
tem and reservoir. After one photon is scattered, the full 
density matrix goes from pf r to due to the trans- 

fer of one atom from the initial modes into the quasi- 
continuum. As the modes are well-separated, the scat- 
tered photon carries sufficient phase information into the 
environment to determine whether the atom recoils from 
the L or the R mode. Assuming that atoms from L and 
R modes scatter into different manifolds of recoil modes, 
the reduced density operator of the two-mode system, 
obtained by tracing out the quasi-continuum modes, be- 
comes 



JV-l 



N t 



(50) 



with Af the normalization factor. In other words the 
system collapses onto an incoherence mixture of the state 
with one atom lost from mode L and the state with one 
atom lost from mode R. 

As an example, we now consider a N particle cat state 

of 



(51) 



and to distinguish cat states with different total atom 
numbers, here we use a new notation 



|* c > = -J=(|iV/4.3iV/4) 



|3JV/4,JV/4», 



(52) 



where \ni,,nji) stands for a Fock basis with atoms in 
left mode and tir in the right. This cat state is collapsed 
after one photon-detection into a statistical mixture of 
two N — 1 particle states, 



„N- 



1 



;*ii + i* 2 )(*2i) 



(53) 



where 
l*i> : 

l*2> 



1 

l 

71 



(|A/4 - 1, 3JV/4) + |37V/4 - 1, N/4)), (54) 
(|JV/4, 3iV/4 - 1) + \3N/4,N/4- 1>). (55) 



The important point here is that each of these states is 
still a good cat-like state, so that while we don't know 
which one the system has collapsed into, we do know that 
the system remains in a good cat-like state. The means 
that the effect of scattering a photon and losing an atom 
may not have a significant detrimental effect on cat-like 
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FIG. 20: The evolution of the projector P ca t under decoher- 
ence die to inelastic atom loss. The curves correspond to the 
dynamically-created cat state (Blue solid line), the perfect 
cat state (Red dashed line) and less extreme cat state IT>21 
(Green dashed-dotted line). The figure plots the estimated 
probability to remain in a cat-like state versus time. Starting 
from 100 atoms, it is estimated that 63 atoms are lost by the 

A 2 

final time t — r^pj; • 



states. Similarly, after a second photon is scattered, the 
system will collapse into a statistical mixture of four N— 2 
particle cat states, and so on. We note, however, that for 
the 'perfect' cat state [\N0) + \0N)] a single scat- 
tering results in a mixture of the states \(N — 1)0) and 
|0(AT — 1)), neither of which is itself a cat-state. This 
means that while the ideal cat state is so fragile that a 
single scattered photon is sufficient to collapse that cat 
onto all-left or all-right states, the cat-like states we are 
considering should be significantly more robust, giving 
them a significant advantage provided they are still suit- 
able for whatever application is desired. 

To show the dynamics under decoherence induced by 
these inelastic condensate losses, in Fig. [201 we plot the 
time evolution of the expectation value of the projector 
P ca t, defined in (|10fl . Three initial states are shown, the 
extreme cat state [\N, 0) + \0N)] /V2, the less extreme 
cat state, [|f + |^, f )] /%/2 and the wave-packet 
cat-states produced via dynamical evolution as in Section 

mm 

Figure |2*U1 is obtained by assuming a single photon is 
detected in each time interval At = rnpfwro ' wriere 
is the number of remaining atoms at time t. After each 
scattering a new density matrix is obtained by applying 
Eq. H50|) . We note that the use of (|5U|I automatically 
weighs the probabilities that the scattering occurs from 
the L or R mode due to the action of the annihilation 
operators being proportional to the square root of the 
mode occupation. The resulting density matrix is then 
diagonalized and written 



(56) 
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FIG. 21: Magnitude of the resulting density matrix elements 
corresponding to Fig. I^OJat time t = 0.3 Jktf ■ 
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FIG. 22: Probability distributions P(n) of the six leading 
eigenmodes of the density matrix shown in Fig. 1211 , 



probability to remain in a cat-like state even after many 
scatterings, and, the dynamically created cat-states are 
similarly robust against decohercncc. We note that in 
the figure out of an initial population of N = 100 atoms, 

A 2 

N = 37 remain at the end-time t = r£rp ■ 

In Fig. [21 we plot the magnitude of the density ma- 
trix after a time of thoid — 0-3 r^pj as in Fig. (|2Tj)l (Blue 
solid line). In this figure, an initial cat state is formed 
at time t ca t via dynamic evolution with no decoherence. 
The state is then held with g — and r = for time 
thoid with photon-scattering acting as a gradual source 
of decoherence. The figure appears to show a collapse of 
the reduced density matrix onto a statistical mixture of 
left-centered and right-centered peaks. However, diago- 
nalization of this matrix reveals that this is not the case. 
Based on numerical diagonalization, we plot the major 
eigenmodes \*f?i) of this density matrix in Fig. [221 From 
the picture, some of these resulting eigenmodes, such as 
l^^'s with i — 1,2,4,5 are still cat-like, in a sense that 
separation between two peaks are larger than the width 
of either peak. While for some others, like IvE^) and 
|\?e)) the separation between peaks are not larger than 
the width of each peak, and thus they are not cat-like. 
As more atoms are scattered, the number of remaining 
atoms decreases, and hence the maximum separation of 
cat-like states will become will necessarily decrease. Be- 
cause we have chosen projector P cat to be fixed- width 
filter, those cat-like state with small separations between 
peaks may not be picked by P ca t- As a result, P cat de- 
creases to almost zero before all atoms are lost, as seen in 
Fig. (|20|l . Lastly, we note that In this section we did not 
show the effect of inelastic decoherence on the tunneling 
dynamics during the formation of the cat state. This 
can be done by the standard quantum jump method in 
Mento-Carlo simulation, which we plan to carry out in 
future work. 



VI. DEPHASING INDUCED BY ATOMIC 
COLLISION 



where denotes the ith eigenmode, with its statis- 
tical weights over the mixture state. The resulting eigen- 
modes |v?»)'s are then analyzed to determine whether or 
not they are sufficiently cat-like. This is accomplished 
by first assigning a value of zero to any state where 
J2n>o\ c n\ 2 S«<ol c «l 2 ' where c n is again the prob- 
ability amplitudes for the state \n). The states which 
survive this test are then assigned a value equal to the 
expectation value of the projector P ca t- The assigned 
values for each eigenmode of the density matrix at time 
t are then averaged together and the results are plotted 
versus time. 

From the figure, we verify that the perfect cat state is 
completely collapsed immediately after one photon scat- 
tering. In contrast, the less extreme cat states has a high 



In this section we consider the effects of atom- atom 
collisions on a cat-state which is 'frozen' by reducing the 
tunneling strength, r, to zero. If during this hold state 
the atomic interaction strength, g, is nonzero, the system 
will undergo a dephasing process j2f| > whose dynam- 
ics is determined by Hamiltonian (|13H . The main effect 
of this dephasing process is that, while it will not collapse 
the coherence of the cat-state, once the state is dephased, 
it cannot evolve back to the initial coherent state. In this 
way, a detection scheme based on revivals will incorrectly 
create the appearance of collape. Therefore, in general, 
the dephasing effect needs to be suppressed by making 
the atomic interaction strength g, and/or the hold time 
small. 

In studying the dephasing effect, first, we note that at 
each period of ir/g, the relative phase of each c„ is ex- 
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FIG. 23: Dephasing process of cat state with nonzero atomic 
interaction. The collapse is measured by P co i — tr{p c p(t)}, 
where p c and p(t) are the densities of the initial cat state and 
dephased state after a hold time thoid respectively. Parame- 
ters are given the same as in Fig. [3] while thoid has units of 



actly restored and thus the many-body states is revived. 
While between two consequent revivals, nonzero relative 
phases are introduced to Fock basis, and the cat state 
will collapse. In Fig. 1231 we show such a rapid dephas- 
ing process, where the cat-state decays with a collapse 
time around t co i — 0.02 / 1 <? | , determined from numerical 
simulations. 

To obtain an analytic estimate of the collapse time 
scale, we assume that our system is prepared in a 
Gaussian-like cat state (©, which is a good approxima- 
tion of the realistic cat state in our system. And then, 
we associate the collapse and revival of our bimodal sys- 
tem with the evolution of one-body two-mode correlation 
function A, 



A = (cat | a^aft | cat 



(57) 



At short time scale, the correlation function l|57|l is ob- 
tained, under the condition of that the Gaussians are well 
peaked around n = ±nrj, to a good approximation [47^ . 



A(i) = A(0)e-( 5CTt > 



(58) 



And therefore, the collapse time t co i can be estimated by 



tr 



1 

2^' 



(59) 



As discuss above, the collisional dephasing will affect 
the revival process. Generally speaking, if the holding 
time thoid is less than collapse time t co i, the system can 
still evolve back to the initial state as shown in Fig. [21 
(a). But under the condition of thoid > tcoi, the cat state 
is dephased, and the coherent state can not be restored. 
Rather it will evolve to a disordered state as shown Fig. 
E|(b). Therefore, for the success of the detection scheme, 
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FIG. 24: Dynamic evolution of dephased cat state after tuning 
on the tunneling. The parameters are chosen as in Fig. 
with a collapse time of t co i = 1.5/r. In the first picture (a), 
the hold time is set to be thoid — tcoi, while for (b), we set 
thoid = 5t co i. It is shown that for (a), the system still revives 
to the initial state, while in (b), the initial state can not be 
restored. 



one should make sure that the dephasing collapse time is 
longer than the hold time and the revival process is only 
affected by decoherence. 



VII. CONCLUSION 

We have considered the double-well BEC system and 
proposed that cat-like states can be formed by switch- 
ing the sign of the interaction strength from positive to 
negative using a Feshbach resonance. We examined the 
possibility of either adiabatic or sudden switching with 
the following results. In the adiabatic case, the ground 
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state evolves into a macroscopic superposition state with 
all of the atoms collectively occupying either the left or 
right well. However, by examining the degeneracy of the 
ground state, we conclude that the cat state is unstable 
against small perturbations and may readily evolve into 
a localized state. For the case of a sudden change in the 
scattering length, we find that an initial condensate state 
evolves dynamically into a double-peaked superposition 
of number states, with one peak corresponding to the 
majority of atoms being in one well, and the other peak 
corresponding to the majority being in the opposite well. 

We have shown that this process is stable against per- 
turbations in the initial state, as well as in the control 
parameters. For this reason we believe that the dynam- 
ical evolution scheme is more likely to be demonstrated 
experimentally. In addition, we have demonstrated that 
continuing the evolution after the cat-state formation re- 
sults in a nearly complete revival of the initial state. 
However, this revival does not occur if the cat-state has 
collapsed into an incoherent mixed state, so that this re- 
vival can be taken as proof of coherence of the initial 
cat-state. We note that the states formed via adiabatic 
switching are maximum 'all or nothing' cat states, while 
those formed dynamically are less distinct cat-like states 
corresponding to a double-peaked distribution in Fock- 
space. 

The effects of decoherence due to spontaneous scat- 
tering of laser photons have also been studied in detail. 
We have determined that decoherence is strongest for 
the case where the single-photon momentum recoil is not 
sufficient to remove the atom from its initial state, due 
to Bose-stimulation as the atom remains in an strongly- 
occupied atomic-field mode. In the opposite regime, 
where the scattering rate is smaller by a factor of N, 
N being the atom number, we find that the maximally 
entangled 'all or nothing' cat states are destroyed by 
the scattering of a single photon, while the dynamically- 
formed less-extreme cat-like states can survive multi- 
ple scatterings without leaving the subspace of cat-like 
states, in analogy to the effect of atomic collision losses 
studied in [2y|. Because of this enhanced survivability, 
the less-extreme cat-states may be more useful for appli- 
cations, providing they are suitable for the desired pro- 
tocol. 

In this context of 'freezing' the cat states for a pro- 
longed time, we have also considered the de-phasing ef- 
fects of atom-atom interactions during this stage, and 
conclude that collisions can mimic the effects of deco- 
herence with regard to a revival- type measurement, how- 
ever with regard to applications, they may or may not be 
detrimental, depending on the type of interaction and/or 
measurements involved. 

We note that in the present proposal, there are three 
different time scales which we have established: the time 
need to create cat state t ca t, the coherence time t co h and 
the collisional dephasing time t co /. First, to generate the 
cat state, the coherence time scale t co h must longer than 
teat- While later, in order to hold the cat state for later 



use in an application, a hold time thoid shorter than t co h is 
required to maintain the macroscopic superposition state. 
And finally, another condition t co i > thoid should be sat- 
isfied to avoid the dephasing quasi-collapse, if such de- 
phasing is detrimental to the desired application. The 
coherence time, t co h, is tunable by adjusting the laser in- 
tensity and detuning, while the dephasing time tdephaae 
can be adjusted if the interaction strength g is modified 
during the hold stage. 



VIII. APPENDIX 

To eliminate the excited state in the master equation 
(|23|l . we note that for a far-off-resonant laser field, the 
excitation rate in BEC is extremely low, which enables 
us to expand the density matrix to a good approximation 
as 



p ~ Poo + Poi + Pio + Pn, 



(60) 



where poo and pn stand for the density matrix of none- 
excitation and one excitation state respectively, while poi 
and pio describe coherence between these two manifolds. 
In fact, it is convenient to define a projection operator 
Pj, which is a projection into a subspace with j atomic 
excitations of the BEC system, such that 



Pi j = I'.l'l'j- 
With this expansion, one obtains 



(61) 



Poo 



= -i /"dVO(r)$t(r)# e (r)p i + 

J d 3 rd 3 r L{r - r')*J(r)*e(r)pii^( r ')*s( r ') 



-if. a, 



Pio 



(62) 



-i j d 3 rn*(r)^l(r)^ g (r)p 00 

- J d 3 rd 3 r L(r - i>)ty(r)$ g (r)fy (r')* e (r')pio 
-iApio, (63) 



Poi — Pioj 



(64) 



and 
Pio 



-i J d 3 rn*(r)¥ e (r)y a (r)p 01 

- J d 3 rd 3 r L(r - r')*£(r)# 9 (r)*J ( r ')# e ( r ') pn 



+H.c. 



(65) 
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In the assumption of low excitation rate, one can adia- ■ f ( f r ^3 r ' L(r — r ) ^t( r ^| r ( r )\pt ( r ')vj/ ( r ')p 10 

batically eliminate poi arL d pio by demanding p 01 ~ J A e 9 9 e 

and pio ~ respectively, yielding (66) 



+ d 3 r—^-pu^l(r)^ g (r) one obtains 



Solving this equation perturbatively up to the order -i , 



Pio = - / d 3 r^*t (r) ^ (r)AM) + + / d 3 r ^M pn ^t (r) ^ (r) 



A " ' ^ J A 

-i J d 3 rd 3 r'd 3 r" M r ~ *^*( r ) ftt( r )$ g ( r )ftt (rO#e(rO^(r")* g (r")p o 

+l /" d 3 rdVd 3 r" ~ O^* ( r ) $t( r )fr g ( r )ftt (^(rO^r^r")^"), (67) 



and poi = Poi" Inserting this result for poiiPio back into 
Eq. (ESI, we find 



Pii = i/ rf 3 rd 3 r "* (r ^ (r,) ^( r) ^ g(r)poo ^t (r0 ^ e(r0 _i J W ;f^^t (r) i s(r) $t (r ' ) i e(r>11 

+ J d 3 rd 3 r'd 3 r"d 3 r" L * {T " ~ ^ ( r M r ) #t (r) ^ g(r)poo ^t ( r "^ e (r"OttI(r')fr g (r')tt g (r")* e (r") 

- /" d 3 rd 3 /d 3 r"d 3 r'" L * (r,/ ~ r/ ^ (r) " (r } ^(r)^ g (r)^t(r w )^ e (r w )pii^:(r )# g (r )* g (r")tt e (r") 

+i y d 3 rdV L(r - r')*£(r)$ fl (r)#t ( r ')$ e ( r ') pu + #. c . ( 68 ) 

I 

Nextt, we apply a similar adiabatic elimination process we normal order the ground state operators in above 
eliminate the excited state matrix p\\. To do that, first equation, and then by setting pn w 0, obtains 



Pii = ij d 3 rd 3 r g!£^ $t (r) $ g(r)poo gt ( r ')* e (r') - i J d 3 rd 3 r *t(r)& g (r')* g (r)* e (r')pn 

+ y d\d?rd\" ^ (r// ~ 2 r 2^ (r)0(r,) ^U^)^ g (^)poo^ g (^O^ g (^O^ g (^")^(r , ') 

-y d 3 rdVd 3 r"dV" ^^lz^M^^^t (r) ^ (r) ^t ( r '")* e (r'")pii*i(r')# g (r')l-t (r ")^ e ( r ") 
+ff.c, (69) 



Again, we solve this equation perturbatively up to order 
■jy, and arrive at 
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Pn = J d 3 rd 3 r'd 3 r" L * {r " ~ 2 r 2^ (r) " (r ) ^(r)# g (r) PO Q^(rO^(rO^(r ,, )# e (r") 

-J d*rd 3 r'dVdV'd 3 r"" ^I^^ x 

^(r'')* g (r")^r>"')$ e (r"") _ J d 3 rd 3 r'dVd 3 r'"d 3 r"" L(r " ^ L *(^^)^^'Mr ) x 
*t(r)*t (r /)tf g(r )# e(r '^ (70 ) 



This express for pn can be simplified by normal order- 
ing. Finally, we insert the normal ordered poiiPoiiPn 
back into Eq. (|62|l . and keep only terms up to i. Mean- 



while, since the number density of atoms in BEC is small, 
we drop any 3-body and above terms. With these two as- 
sumptions, we obtain the resulting master equation 12711 . 



Ps (t) = - iHsPs (t)-^r J d 3 r¥ g (r)i> g (r)p s (t)- 

I Mrl O(r)£(r-r0n*(r0 ^t (r) ^t (r) ^ g(r0 ^ (r0/3s(t) + 
j d 3 rd 3 r ^*(r)X(^-rOO(rO # t (r) ^ (r)p3(f) ^t( r0 ^ g(r0+g ,, 



where we have written p s = poo- 
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